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ABSTRACT 

To maximize the number of planet detections by increasing efficiency, current microlensing follow-up obser¬ 
vation experiments are focusing on high-magnification events to search for planet-induced perturbations near 
the peak of lensing light curves. It was known that by monitoring high-magnification events, it is possible to 
detect multiplicity signatures of planetary systems. However, it was believed that the interpretation of the sig¬ 
nals and the characterization of the detected multiple-planet systems would be difficult due to the complexity of 
the magnification pattern in the central region combined with the large number of lensing parameters required 
to model multiple-planet systems. In this paper, we demonstrate that in many cases the central planetary per¬ 
turbations induced by multiple planets can be well approximated by the superposition of the single planetary 
perturbations where the individual planet-primary pairs act as independent binary lens systems (binary super¬ 
position). The validity of the binary-superposition approximation implies that the analysis of perturbations 
induced by multiple planets can be greatly simplified because the anomalies produced by the individual planet 
components can be investigated separately by using relatively much simpler single-planetary analysis, and thus 
enables better characterization of these systems. 

Subject headings: planetary systems - planets and satellites; general - gravitatinal lensing 


1. INTRODUCTION 

Microlensing is one of the most important methods that can 
detect and characterize extrasolar planets (see the review of 
Perryman 2000). Microlensing planet detection is possible 
because planets can induce perturbations to the standard lens¬ 
ing light curves produced by primary stars. Once the per¬ 
turbation is detected and analyzed, one can determine the 
mass ratio, q, and the projected separation, s (normalized by 
the Einstein ring radius r p), between the planet and host star 
dMao & Paczvnsk3ll991h . Recently, a cle ar-cut microlensin g 
detection of an exoplanet was reported bv iBond et alJ (120041) . 

Planet detection via microlensing is observationally chal¬ 
lenging. One of the most important difficulties in detecting 
planets via microlensing lies in the fact that planet-induced 
perturbations last for a short period of time. Eor a Jupiter- 
mass planet, the duration is only a few days and it decreases 
as oc y/q. To achieve high monitoring frequency required for 
planet detections, current lensing experiments are employ¬ 
ing early warning system to issue alerts of ongoing events 
in the early stage of le nsing magnification dSoszvnski et alJ 
1200 It iBond et ^120011) and follow-up obs ervation programs 
to intensively monitor the alerted e vents dBond et al .1120021 
iPark et al1l2004t ICassan et alJl2004]) . However, follow-up is 
generally done with small field-of-view instrument, and thus 
events should be monitored sequentially. As a result, only a 
handful number of events can be followed at any given time, 
limiting the number of planet detections. 

An observational strategy that can dramatically in- 
crease the planet detect ion efficiency was proposed by 
iGriest & Safizadell (119981) . When a microlensing event is 
caused by a star possessing a planet, two sets of caustics are 
produced. Among them, one is located away from the pri¬ 
mary lens (planetary caustic) and the other is located close to 
the primary lens (central caustic). The location of the plan¬ 
etary caustic varies depending on the planetary separation s, 
which is not known, and thus it is impossible to predict the 
time of planetary perturbation in advance. On the other hand, 
the central caustic is always located very close to the primary 


lens, and thus the central perturbation occurs near the peak of 
high-magnification events. Therefore, by focusing on high- 
magnification events, it is possible to dramatically increase 
the planet detection efficiency, enabling one to maximize the 
number of planet detections with a limited use of resources 
and time dHan & Kiml200 lURattenburv et alJ2002HYoo et alJ 

Hool. 

An ad ditional use of high magnificat ion events was no¬ 
ticed by IlGaudi, Naber, & S^ketll ill998ft . They pointed out 
that multiple planets with separations ~ 0.6 - 1.6 of the Ein¬ 
stein ring radius significantly affect the central region of 
the magnification pattern regardless of the orientation and 
thus microlensing can be used to detect multiple-planet sys¬ 
tems. They noted, however, that characterizing the detected 
multiple-planet systems by analyzing the central perturbations 
would be difficult due to the complexity of the magnification 
pattern combined with the large number of lensing parameters 
required to model multiple-planet systems. 

In this paper, we demonstrate that in many cases the cen¬ 
tral planetary perturbations induced by multiple planets can 
be well approximated by the superposition of the single¬ 
planetary perturbations where the individual planet-primary 
parrs act as independent binary lens systems (binary superpo¬ 
sition). The validity of the binary-superposition approxima¬ 
tion implies that a simple single-planet lensing model is pos¬ 
sible for the description of the anomalies produced by the in¬ 
dividual planet components, enabling better characterization 
of these systems. 

The layout of the paper is as follows. In § 2, we describe the 
multiple-planetary lensing and magnification pattern in the 
central region. In § 3, we illustrate the validity of th binary- 
superposition approximation in describing the central magni¬ 
fication pattern of multiple-planet systems. In § 4, we discuss 
the usefulness of the superposition approximation in the inter¬ 
pretation of the multiple planetary signals. 


2. MULTIPLE PLANETARY LENSING 
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Fig. 1.— Contour maps of magnification excess from the single-mass leasing as a function of the source position (^. 77 ) for a lens system containing two 
planets. The parameters of planet 1 (located on the left-side ^ axis) are held fixed at q[ = 0.001, = 1.4, while the projected separation S 2 and the angle between 

the axes, A0, are varied for a second planet with ^2 = 0.0003. Greyscale is used to represent positive (bright) and negative (dark) deviation regions. The two sets 
of contours drawn in black and white are for the maps constructed by exact triple leasing and binary-superposition approximation, respectively. Contours are 
drawn at the levels of |e| = 2%, 5%, and 10%. The straight line with an an'ow in each panel is the source trajectory where the resulting light curve is presented 
in the corresponding panel of Fig.|3| Lengths are scaled by the Einstein ring radius 6^. The coordinates are centered at the effective position of the primary lens 
[see eq. () for definition]. The orientation is such that /S.9 = 0° implies the two planets are on the same side. 


Fig. 2.— Enlargement of the excess magnification maps presented in Fig.Q The figures drawn in solid black and white lines are the caustics for the cases of 
the exact triple lensing and binary-superposition approximation, respectively. Contours and greyscale are drawn by the same way as in Fig.lTI 


The equation of lens mapping from the lens plane to the 
source plane (lens equation) of an N point masses is expressed 
as 
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where C = ^ + = xlj + iyLj, and z = x + iy are the com¬ 

plex notations of the source, lens, and image positions, re¬ 
spectively, z denotes the complex conjugate of z, nij are the 
masses of the individual lens components, M = '^^e 

total mass of the system, and thus irij/M represent the mass 
fractions of the individual lens components. Here all angles 
are normalized to the Einstein ring radius of the total mass of 
the system, M, i.e. 
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where Dql and Dqs are the distances to the lens and source, 
respectively. The lensing process conserves the source surface 
brightness, and thus the magnifications A,- of the individual 
images correspond to the ratios between the areas of the im¬ 
ages and source. Lor an infinitesimally small source element, 
the magnification is. 
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( 3 ) 


The total magnification is the sum over all images, A = X); A,. 

Lor a single lens (N = 1), the lens equation is simply in¬ 
verted to solve the image positions (.r,y) and magnifications 
for given lens {xi^yi) and source (^, 77 ) positions. This yields 
the familiar result that there are two images with magnifi¬ 
cations and separations from the lens of A± = 0.5(A ± 1) 
and ui_± = \z±-Zl\ = 0.5[m±(m^-I-4)*/^], respectively, where 
M = |C“Zl| is the separation between the lens and source and 
A = A+-l-A_ = (m^-I-2)/[m(m^-I-4)*/^] is the total magnification. 

A planetary lensing with a single planet is described by the 
formalism of a binary {N = 2) lens. In this case, the lens equa¬ 
tion cannot be inverted algebraically. However, it can be ex¬ 
pressed as a fifth-order polynomial in z and the image posi¬ 
tions are then obtained by numerically solving the polyno¬ 
mial (iWitt & Maoll 19951) . One important characteristic of bi¬ 
nary lensing is the formation of caustics, which represent the 
set of source positions at which the magnification of a point 
source becomes infinite. Planetary perturbations on lensing 
light curves occur when the source approaches close to the 
caustics. The location and size of these caustics depend on 
the projected separation s and the mass ratio q. Lor a plane¬ 
tary case, there exist two sets of disconnected caustics. The 
planetary caustic(s) is (are) located away from the primary 
star on or very close to the star-planet axis with a separation 


from the primary star of ^pc ~ s - 1/s. The central caustic, 
on the other hand, is located c lose to the primary lens with 
a size of A^cc — 4^/(s- 1/s)^ (IChung et al.ll2005li . Caustics 
are located within the Einstein ring when the planetary sep¬ 
aration is within the range of 0.6 < s < 1.6. The size of the 
caustic, which is directly proportional to the planet detection 
efficiency, is maximized when the planet is located within 
this range, and thus t his range is referred as “lensing zone” 
jr.ould&Loehl1992h . 

Lor a multiple-lens system, the lens equation is equiva¬ 
lent to a polynomial with an order of {N^+ 1). Therefore, a 
multiple-lens system produces a maximum of -t 1 and a 
minimum of A-I-1 images, with the number of images c hang¬ 
ing by multiple of two as the source crosses a caustic dRhiel 
fT9^ . Lor a system with two planets (and thus A = 3), there 
are thus a maximum of 10 images and a minimum of 4 im¬ 
ages. Unlike the caustics of binary lensing, those of multiple 
lensing can exhibit self-intersecting and nesting. 

3. BINARY-SUPERPOSITION APPROXIMATION 

The mass ratio of a planet to its primary star is very small, 
and thus the lensing behavi or of planet-induced anomalies can 
be treated as perturbati on llGould & LoebllT9^lBozzall 19991 
lAsadall200^ lDominildl1999t lAnll2005h . Due to the pertur¬ 
bative nature of planetary anomalies, it was known that the 
magnification pattern in the region around planetary caustics 
of multiple-planet systems can be described by the superposi¬ 
tion of those of the single-planet systems where the individual 
planet-primary p airs acts as independent binary lens systems 
(iHan et d EOOI!), i.e.. 


A^-1 

Amul ^ Abs = ^ ^ Abi,/~Ao, (4) 

(=1 

where Amui, Abs, Abi,;, and Aq represent the magnifications of 
the exact multiple lensing, binary-superposition approxima¬ 
tion, binary lensing with the pairs of the primary and indi¬ 
vidual planets, and single-mass lensing of the primary alone, 
respectively. By contrast, it was believed that the binary- 
superposition approximation would not be adequate to de¬ 
scribe the magnification pattern in the central region because 
the central caustics produced by the individual planet com¬ 
ponents reside at the same central region and thus non-linear 
interference of the perturbations would be large. 

Unlike this belief about the magnification pattern in the cen¬ 
tral perturbation region, we find that non-linear interference 
between the perturbations produced by the individual planets 
of a multiple-planet system is important only in a small con¬ 
fined region very close to the central caustics. This implies 
that in many cases binary-superposition approximation can 
also be used for the description of the magnification pattern 
in the central perturbation region. 
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Fig. 3.— Greyscale maps of the difference between the magnification excesses of the triple lensing and binary approximation, i.e. Ac = Cbs^ctri = (^bs“ 
^tri)Mo- Greyscales brighter and darker than background represent the region where Atn > Abs and vice versa, respectively, and the scale changes at |ebsl = 2% 
and 5%. 


Fig. 4.— Example microlensing light curves produced by multiple-planet systems. The parameters of the planets are identical to those in Fig.Q The source 
trajectories responsible for the events are marked also in Fig.0 For all events, the impact parameter of the trajectories is uq = 0.03. In each panel, there are five 
curves, where the curves drawn in thick solid black, thin magenta, cyan, and red lines are those of exact triple lensing, binary lensing with the pairs of the primary 
and the individual planets, and binary-superposition approximation, respectively, and the dashed curve is that of a single-mass lensing of the primary alone. The 
curve in the lower part of each panel shows the difference between the excesses of the exact triple lensing and binary-superposition approximation. 


To demonstrate the validity of the binary-superposition ap¬ 
proximation in the central region, we construct a set of magni¬ 
fication excess maps of example multiple-planet systems con¬ 
taining two planets with various orientations. The magnifica¬ 
tion excess represents the deviation of the magnihcation from 
the single-mass lensing as a function of the source position 
(^, rj), and it is computed by 


Aai-Ao 
Ctrl - —T- 

Ao 


( 5 ) 


where Atd is the magnification of the triple (primary star plus 
two planets) lensing. 

In Figure n we present the constructed contour (drawn by 
black lines) maps of magnification excess. In the map, the 
parameters of planet 1 are held hxed at = 0.001 and si = 1.4, 
while the projected separation S 2 and the angle between the 
position vectors to the individual planets from the primary star 
(orientation angle), A9, are varied for a second planet with 
q 2 = 0.0003. Greyscale is used to represent positive (bright, 
Etri > 0) and negative (dark) deviation regions. Also drawn are 
the contours (drawn in white lines) of magnification excess 


based on binary-superposition approximation, i.e.. 


Ebs = 


E 


Abi,/ Aq 
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For the maps based on binary superposition, we consider 
slight shift of the effective lensing position of the primary star 

i Zf <,) toward the individu al planets (zip.) with an amount of 
Di Stefano & Maoll996h 
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To better show the magnification pattern in the very central 
region and the detailed caustic structure, we enlarge the maps 
and presented them in Figure |2l In each map, the figures 
drawn in thick black and white lines represent the caustics for 
the cases of the exact triple lensing and binary superposition, 
respectively. 

From the comparison of the excess maps constructed by the 
exact triple lensing and binary-superposition approximation, 
one hnds that binary superposition is a good approximation 
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in most of the central perturbation region as demonstrated by 
the good match between the two sets of contours. Slight devi¬ 
ation of the binary-superposition approximation from the ex¬ 
act lensing magnification occurs (a) in a small region very 
close to the central caustics and (b) in the narrow regions 
along the primary-planet axes. This can be better seen in 
Ligure|3 where we present the greyscale maps for the dif¬ 
ference between the excesses of the triple lensing and binary- 
superposition approximation, i.e., 

= Ctrl ~ Cbs- (8) 

The difference in the region close to the central caustics is 
caused by the non-linear interference between the perturba¬ 
tions produced by the individual planets. On the other hand, 
the difference along the primary-planet axes is caused by the 
slight positional shift of the triple-lensing caustics due to the 
introduction of an additional planet. However, the area of the 
deviation region, in general, is much smaller than the total 
area of the perturbation region, and thus the binary superposi¬ 
tion approximation is able to well describe most part of plan¬ 
etary anomalies in lensing light curves. This can be seen in 
Ligure|4] where we present example light curves of multiple¬ 
planetary lensing events and compare them to those obtained 
by binary superposition. 



^e (deg) 


Fig. 5.— Dependence of the binary-supeiposition validity (7^) on the 
angle between the msition vectors to the component planets from the host 
star (A6). See eq. for the definition of 7^.. To see the dependence only on 
A6, we fix (^ 1 , ) and (^21 •* 2 ) imd their values are marked in the panel. 

The lensing behavior of a multiple planetary system is de¬ 
termined by many parameters including q\, q 2 , si, S 2 , and 
A6, and binary-superposition might not be a good approxi¬ 
mation in some space region of these parameters. We there¬ 
fore investigate the region of parameter space where binary- 
superposition is a poor approximation. Due to the numerous¬ 
ness of the parameters, we choose a method of investigation 
where we inquire the validity of the approximation on the in¬ 
dividual parameters by varying one parameter and fixing other 
parameters. The validity of the approximation is quantified by 
the ratio of 

7^=^, (9) 
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Fig. 6.— A contour map of the binary-superposition validity (TZ) as a 
function of the planet mass ratios of a multiple planetary system composed 
of two planets. To see the dependence only on the mass ratios, we fix ii, S 2 
and AO and their values are marked in the panel. 


where and Eac represent the area of the central pertur¬ 
bation region and the area of the difference region between 
the exact triple lensing and the binary-superposition approxi¬ 
mation, respectively. In some cases, the perturbation regions 
caused by the planetary and central caustics are connected, 
making the boundary between the two regions ambiguous. 
We thus define the central perturbation region as the region 
within the lens-source impact parameter of m = 0.1 (corre¬ 
sponding to the region with magnifications A > 10). The areas 
Eae and E^ are computed by setting the threshold values of 
\Ae,h \ = 5% an^er/il = 5%, respectively. 

In Ligures|3|3 and0 we present the dependence of TZ on 
the orientation angle ( AO), the mass ratios of the component 
planets (qi and q 2 ), and the separations to them (ii and S 2 ), re¬ 
spectively. Lrom the variation of TZ depending on the orienta¬ 
tion angle, we find that the difference between the exact lens¬ 
ing and binary-superposition approximation becomes bigger 
as AO decreases and thus the two planets are located closer 
to each other. We interpret this tendency as the increase of 
the non-linear interference between the perturbations caused 
by the two planets as the separation between them decreases. 
Lrom the dependence on the mass ratio and separation, we 
find that binary superposition becomes a poor approxima¬ 
tion as either the planet mass increases or the separation ap¬ 
proaches to unity. These tendencies are the natural results of 
the breakdown of the perturbation treatment for companions 
with high mass ratios and (or) separations of s ~ 1, because 
the perturbation treatment is valid when q-^1 and |i-11 ^ q'. 
Besides these extreme regions of parameter space, however, 
we find that R < 5% in most regions, implying that binary- 
superposition approximation well describes the lensing be¬ 
havior of most multiple planetary systems. 

4. IMPLICATION AND CONCLUSION 

In the previous section, we demonstrated that lensing mag¬ 
nification patterns of multiple-planet systems can be de¬ 
scribed by using binary-superposition approximation not only 
in the region around planetary caustics but also in most part of 
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Fig. 7.— A contour map of the binary-superposition validity (TZ) as a 
function of the separations to the planets from the host star of a multiple 
planetary system composed of two planets. To see the dependence only on 
the separations, we fix qi, q 2 and AO and their values are marked in the panel. 


the central perturbation region. In this section, we discuss the 
importance of the binary-superposition approximation in the 
analysis and characterization of multiple-planet systems to be 
detected via microlensing. 

Exact description of lensing behavior of events caused by 
multiple-planet systems requires a large number of parame¬ 
ters. Even for the simplest case of a two-planet system, the 
number of parameters is ten, including four single-lensing 
parameters of the Einstein timescale fg, time of the closest 
lens-source approach to, lens-source impact parameter uq, and 
blended light fraction fs, and another four parameters of the 
two sets of the planetary separation and mass ratio, 
and {s 2 ,qi), and the source trajectory angle a and the orien¬ 
tation angle A0. As a result, it was believed that analyzing 


anomalies produced by multiple planets would be a daunting 
task. 

With the validity of binary-superposition approximation, 
however, the analysis can be greatly simplified because the 
anomalies induced by the individual planets can be inves¬ 
tigated separately by using relatively much simpler single¬ 
planetary lensing analysis. Anomalies for which this type 
of analysis is directly applicable are those where the per¬ 
turbations induced by the individual planets are well sepa¬ 
rated, e.g., the anomalies in the lensing light curves presented 
in Eigureglwith fe, A6») = (0.7,60°), (0.7,120°), (0.9,60°), 
(0.9,120°), (1.3,60°), and (1.3,120°). 

In some cases, the region of perturbations caused by the 
individual planets are located close together and thus part of 
the anomalies in lensing light curves can be blended together. 
However, even in these cases, the non-linear interference be¬ 
tween the anomalies is important only in small confined re¬ 
gions and thus the superposition approximation would still be 
valid for a large portion of the anomalies, allowing rough es¬ 
timation of the separations and mass ratios of the individual 
planets. Once these rough parameters are determined, then 
fine tuning of the parameters by using exact multiple-planet 
analysis will be possible by exploring the parameter space that 
was greatly narrowed down. 

However, care is required for the case when the two per¬ 
turbations caused by the individual planets happen to locate 
at the same position (or very close to each other). This case 
occurs when the two planets are aligned (A6* = 180°) or anti¬ 
aligned (A6* = 0°). In this case, the anomaly appears to 
be caused by a single planet and thus naive analysis of the 
anomalies can result in wrong characterization of the planet 
system. However, this type of anomalies will be very rare. 

We would like to thank A. Gould for making useful com¬ 
ments on the paper. This work was supported by the As- 
trophysical Research Center for the Structure and Evolution 
of the Cosmos (ARCSEC") of Korea Science & Engineer¬ 
ing Eoundation (KOSEE) through Science Research Program 
(SRC) program. 
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